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Abstract 
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We study the half advanced and half retarded Wheeler Green func- 
tion and its relation to Feynman propagators. First for massless equa- 
tion. Then, for Klein-Gordon equations with arbitrary mass parame- 
ters; real, imaginary or complex. In all cases the Wheeler propagator 
lacks an on-shell free propagation. The Wheeler function has support 
inside the light-cone ( whatever the mass ) . The associated vacuum is 
symmetric with respect to annihilation and creation operators. 

We show with some examples that perturbative unitarity holds, 
whatever the mass ( real or complex ). Some possible applications are 
discussed. 

PACS: 10. 14. 14.80-j 14.80.Pb 
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1 Introduction 

More than half a century ago, J. A. Wheeler and R. P. Feynman plubhshed 
a work [1] in which they represented electromagnetic interactions by means 
of a half advanced and half retarded potential. The charged medium was 
supossed to be a perfect absorber, so that no radiation could possibly scape 
the system. 

We are going to call this kind of potential a "Wheeler function" ( or 
propagator ), although it had been used before by P. A. M. Dirac [2] when 
trying to avoid some run-away solutions. Later, in 1949, J. A. Wheeler and 
R. P. Feynman showed that, in spite of the fact that it contains an advanced 
part, the results do no contradict causahty [3]. 

Of course, the success of QED and renormalization theory made soon 
unnecesary or not advisable, to follow that line of research ( at least for 
electromagnetism ). 

One of the distinctive characteristics of the Green function used in refer- 
ences [1, 2, 3] is its lack of asymptotic free waves. This is the reason behind 
the choice of a "perfect absorber" for the medium through which the field 
propagates. As the quantization of free waves is associated to free particles. 
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the above mentioned feature of Wheeler functions imply that no free quan- 
tum of the field can ever be observed. Nevertheless, we are now habituated 
to the existence of confined particles. They do not manifiest themselves as 
free entities. 

We can give some examples ( outside QCD ) where such a behaviour can 
be present. 

A Lorentz- invariant higher order equation can be decompossed into Klein- 
Gordon factors, but the corresponding mass parameters need not be real. For 

instance, the equation: 

(□^ + m^'^cp= (n + im^^ (□ - iw?) ip = Q (1) 
gives rise to a pair of constituent fields [4] obeying: 

(□ ± im^) = (2) 

Any solution of eq.(2) blows- up asymptotically. We can say that the 
corresponding fields should be forbidden to appear asymptotically as free 
waves. Therefore, they should have a Wheeler function as propagator [5]. 

Equations similar to (1), or more general: 



4 



± m^") ip^O (3) 

appear in a natural way in supersymmetric models for higher dimensional 
spaces [6]. 

Another example is provided by fields obeying Klein- Gordon equations 
with the wrong sign of the mass term. A careful analysis shows that the 
propagator should be a Wheeler function [7, 8]. Accordingly no tachyon can 
ever be observed as a free particle. They can only exist as "mediators" of 
interactions. 

To define the propagators in a proper way, we have to solve the equations 
for the Green functions, with suitable boundary conditions. 
For the wave equation: 

DG{x) = S{x) 
a Fourier transformation gives: 

G{p) = {f - pI) ~" ^ ip^p^-' ^ (5) 
Of course, it is necessary to specify the nature of the singularity. Different 
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determinations imply different types of Green functions. For tlie classical 
solution of (4) it is natural to use the retarded function (Grt)- It corresponds 
to the propagation towards the future of the effect produced by the sources. 
This function can be obtained by means of a Fourier transform of (5) in 
which the po integration is taken along a path from — oo to +00 , leaving the 
poles to the right. In practice, we add to po a small positive imaginary part: 



Grt{p)= [p^-{po + lOf 

(p^ - Po - ^0 sgnpo) = {P -iO sgnpoY^ (6) 
The advanced solution is the complex conjugate of (6): 

Gadip) = (p^ - Po + ^0 sgnpt^ = (P + iO sgnpo)'^ (7) 
For the Feynman propagator we have to add a small imaginary part to 
P ( not just to Po ) '■ 

G±{p) = {P±tO)-' (8) 
And, in the massive case: 

Gilp) = (P + m^ ±i0)~^ (9) 
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The Wheeler function is half advanced and half retarded. It is easy to 
see that it is also half Feynman and half its conjugate ( we will not use any 
index for the Wheeler propagator): 

= + (10) 

On the real axis, the Wheeler function coincides with Cauchy's "principal 
value" Green function, which is known to be zero on the mass-shell ( no free 
waves ). 

We can write: 

G±ip) = G{p) ± inS (P + m^) 

where: 

m5{P + m')=\G^{p)-\G_{p) (12) 

Equation (11) is a decomposition of the Feynman function into two terms. 
The first one only contains virtual propagation. The second one is a Lorentz 
invariant solution of the homogeneous equation representing the free particle. 

To perform convolution integrations in p-space, we will utilize the method 
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(11) 



presented in reference [9]. Essencially, it consists in the use of the Bochner 
theorem for the reduction of the Fourier transform to a Hankel transform. 
The nucleus of this transformation is made to correspond to an arbitrary 
number of dimensions v , taken as a free parameter. In this way, starting with 
a given propagator in p-space, we get a function in x-space whose singularity 
at the origin depends analytically on i/ . It exists then a range of values ( of 
V ) such that the product of Green functions is allowed and well determined. 
In x-space we define: 



By means of a "Wick rotation" , the po-integration can be made to run 
along the imaginary axis, without crossing any pole. Mathematically, we 
perform a dilatation: 




The Fourier transform of the massless Feynman function is: 




1 , 



(13) 



A subsequent continuation to a = i produces the transformation: 
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The new quadratic forms are Euclidean and Bochner theorem [10] tells 
that the Fourier transformation reduces to: 



oo 



J-{(P-iO)-n(x) = 4 (14) 



where J^a is a Bessel function of the first kind and order a . 
In eq.(13) we sec that a Wick rotation in p-space (a — )■ i ) implies an 
anti-Wick rotation in x-space {a~^ — > — i ). We must then choose: 

x = (g + io)^ 

Note that the imaginary unit in eq.(14) ( r.h.s. ), is due to the transfor- 
mation dpo — > idpQ 

Prom ref.[ll] we take: 



I.e.: 



More general, for a function f{P ± iO), we obtain: 



(15) 



:F{f{P±tO)}{x) = T^ Jdyy^f{y')j^_,{xy) (16) 



where x = {Q ^iO)^. 

The right hand side of eq.(16) is a Hankel transform of the function /(y^) 
[12]. 

2 Massless case 

a) Fourier transforms 

With the procedures described in section 1, we can obtain the Fourier 
transforms of general massless Feynman functions, defined as (P ± iO)". 
From (18) and ref.[ll], we get ( compare with ref.[13] ): 

T{{P ± ^0)"} (x) = ^i2^«+l£i^±|l (Q ^ zO)-"-^ (17) 
The exponent oi Q ^ iO can be deduced by dimensional considerations. 
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Furthermore, if we interchange the quadratic forms P ■H- Q and write J-" 
for then eq.(17) is still vahd. 

We define the massless Wheeler propagator as: 



( We will not use any index when a Wheeler function is meant ) 
The Fourier transform of (18) is ( cf. eq.(17) ): 



(18) 



r(a + f 

r(-«) 



-(g + ^o)-"-^--(g-^o)- 



(19) 



But we also have the relation ( valid for any quadratic form [13]) 



(g ± i^f = + e±^^^g^ 



(20) 



where 




{-QY g<o 
g > 
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So that we can write (19) in the form: 

J- {P"} {x) = 22"+f li^^illsmTT (a + ^"j Ql""^ (21) 

r (— a) V 2/ 

Equation (21) shows another interesting property of Wheeler functions. 
They are real and have support inside the light-cone of the coordinates. 

Furthermore, for a — —1, the trigonometric function tends to zero for i/ — > 4 

1— - 

, but Q_ ^ has a pole aX v = A with residue 5{Q) [14]. Then: 

Iimj-{P-I}(x) = 5(g) (22) 

In four dimensions the massless Wheeler function is concentrated on the 
hght cone. Prom eq.(21) we obtain: 

22A+ir (A + 1) r (a + p--^-^ (23) 

where the relation 

T{z)T{l -z)^ 

SmTTZ 

has been used. 

From (17) we can also get: 
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•F-^K}(p) = -2^^+^r(A + i)r(A+0 X 



—A—— TT — A— — 

cosnX P+ ^ + cos—u P_ ^ 



Prom relation (12) we find: 



2mS{P) = (P - iO)"^ - (P + ^0) 



And 



J-{5(P)}(x) 



22-2 



27r 



r x-1 [(g + io)'-^ + (g-io)'-^ 



(5+ ^ - cos-vQ_ 2 



Note that eq.(21) allows us to obtain the advanced and retarded 
nents of the Wheeler functions in x-space. 

Grt{x) = 2e(t)^{P"} {x) Gad{x) = 2Q{-t)T{P^} {x) 

where Q{t) is Heaviside's step function, 
b) Convolutions 
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The well-known convolution theorem asserts that the Fourier transform of 
a convolution product is the product of the Fourier transform of each factor: 

f{p) * 9{P) = c7-^ {-^{/(P)} {^)^{9m {x)} (p) 

c=(27r)^ (28) 

Formula (28) can be taken as a definition of the operation *. The product 
of distributions inside the curly brackets, can be taken in a suitable range of 
and analytically extended to other values. 
For Feynman propagators (28) gives: 

[a + A 

(p-ior*(p-tor = ic2-^—^, — ^ X 

{—oi) 



r -2a - ^ 



[P - iO)'"+2 (29) 



r {2a + 1/) 

An analogous equation can be obtained by changing the sign of the imag- 
inary unit in (29). The convolution of two Feynman functions of the same 
kind (-1- or -) gives another Feynman function of the same kind. 

For Wheeler propagators, eqs.(21), (23) and (28) give: 

p. p. + t) r {-2a - f ) 

r2 {-a) T {2a + z/) 

tgn [a + P^"+t (30) 
14 



So that the convolution of two Wheeler functions gives another Wheeler 
function. 

For the wave equation we choose a — —1: 



{P - iOy^ ■ {P - iOy^ = 2ia{u) {P - iO)^'^ (31) 



(P + iOy'^ ■ {P + iOy'^ = -2ia{u) {P + iO)^ 



p-Up-^^a{v) tgn{^-l)P^^-^ 



(32) 



(33) 



where 



V2 J r(u-2) 



(u-2) 

Eqs.(31) and (32) have a pole for i/ ^ 4 (the usual ultraviolet divergence), 
while (33) is well determined in that limit. 

The convolution of 5{P) function can be found with the help of eqs.(23), 

(24), (27) and (28): 



7r^6{P) * 6{P) = a(i/) tgTr (^^ - 1 



^-2 TT i^~2 

PI — COS—fP^ 



(34) 



A comparison with eq.(31) shows that: 



7r25(P) * 5{P) = sgnP ■ P'^ * P'^ 



(35) 
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This relation implies: 



(P - iO)"^ * (P + iO)"^ + n'^S{P) * 5{P) = 2e(P) P"^ * P"^ 

The convolution of two Feynman propagators of different kinds has sup- 
port outside the light-cone in p-space. 

3 Massive case 

A bradyon field obeys a normal Klein-Gordon equation. Its Feynman prop- 
agator is given by eq.(9). The Wheeler function is: 



(P + m^) ^ = ^ (P + + iO) V ^ (P + - io) ^ (36) 
The on-shell 5-function, solution of the homogeneous equation is: 



1 
27ri 



(P + + iO) ^ - (P + - iO) ^ 



(37) 



a) Fourier transforms 

To find the Fourier transform of the Feynman propagators we use eq.(16) 
and ref.[ll] ( p. 687 - 6.566 - 2 ): 



'^1 (x) = Tim^~^ (gTiO)5(^-5);C|_i [m{Q^iO)^ 
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where /Cq, is a Bessel function of the third kind and order a. (For the 
definition of different Bessel functions we follow ref.[ll], p. 951 - 8.40). Note 
also that with the more general formula of ref.[ll] (p. 687 - 6.565 - 4), we could 
find the Fourier transform of arbitrary powers of the Feynman propagator. 

Using the relation (20), we can write: 



|(P + m^±i0) ^|(a;) = =Fim2 ^ 



where /3 = 1 for the upper sign and /3 = 2 for the lower sign. 
It is now easy to find the transforms of (36) and (37): 



(38) 



-1 



TT 



(a;) = -mt-iQ!^ ^'j,_.{mQ 



(39) 



1 ii 
-— m2 

TT 



Q: 



/C^_i imQ 



(40) 



Also for the massive case, the Wheeler function is zero outside the light- 



cone. 



With some sligth changes in notation, we can check that (38) and (40) 
coincide with the results of ref.[13]. 
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b) Convolutions 

We are going to follow the procedure already used in section 2. 

For the convolution of a massive Feynman propagator with a massless 

one, eqs.(17), (38) and (28) give: 

(P + - io) ~\{P- iO)~^ = - l) x 

|(g + iO)^('-^) /C|_i [m {Q + iO)^] I (p) 
And with the help of eq.(16) and ref.[ll] (p.643 - 6.576 - 3), we get: 

{p + m^- ^0)-' * (P - = ^2^^^r - l) T (2 - x 

where F{a, b, c; z) is Gauss hypergeometric function. 

For two Wheeler propagators and two on-shell 5-functions we can follow 
the same method. It is then not difficult to prove the relation (compare with 
(35)): 

n^S (P + m^) * 5 (P) = sgnP ■ (P + m^)"^ * P'^ 
Or, more general: 
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7r^(5 (P + mfj*5(^P + mfj = sgnP ■ (^P + mfj % (P + mfj ^ (42) 



4 Tachyons 

A tachyon field obeys a Klein- Gordon equation with the wrong sign of the 
"mass" term. The Green function is an inverse of P — //^ (we use //^ = — 
for the mass of the tachyon). 

Although it is not easy to see what a Feynman propagator should be when 
the inverse of P — n'^—p^ —p^ — jj? has a pair of imaginary poles {if p"^ < jJ^), 
we may nevertheless define: 

(P - //^ ± iO)~^ = - (-P + A*^ T iO)~^ 
Or, introducing the "dual" quadratic form: 

+p^-p = pI- f 

In other words, instead of propagators with the wrong sign of the mass, 
we can say that we have propagators with the wrong sign of the metric. 
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To find the Fourier transform, we recall (see section 2) that we have an 
"i" for the Wick rotation of po- For the dual metric we have three Wick 
rotations. The three factors dp contribute with — —i. We also note that 

'^Q+ = Q- and ^Q- = Q+. So that (compare with (38)): 



tUp- ii^±i0) [ (x) = ±i//5-i 



Ql 



2 ^Kf-) 



(43) 



where ^ = 2 for the upper sign and ^ = 1 for the lower sign. 
The real part of (43) is: 

{ (P - /.^ + ZO) } {X) + I (P - _ -1 1 



ReJ=\^[P-^i'±^Q) ^{x) = -^i-^-'Ql' (/.g+) (44) 

This real part has support outside the light-cone, while for bradyons, the 
real part of the Feynman propagator is zero for x'^ space- like (cf. eq.(39)). 

We will now show that (44) is not the Wheeler propagator for the tachyon. 
To this aim, we go back to the original definition. Namely, a half retarded 
and half advanced propagator. 
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fx' 



-1 1 



Ad 2 



The Fourier transform is: 



1 1 



Ad 



1 1 



(271)^^/' 



We will first evaluate the advanced part of (46): 



(45) 



(46) 



-,—ipoxo 



-Po- 1^"^ 



(47) 



Where the path of integration runs parallel to the real axis and below 
both poles of the integrand. For xq < the path can be closed on the upper 
half plane of po- The contribution will be the residues at the poles: 



Hi 



Po = i^*^ = =t^ Y IJ'^ — P^ , if P^ < A*^ 



2uj 2uj 
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-27r 



sinujXQ 



r\(^2 2\ , shu Xq / 2 -2\ 

Q[p - /X j H [l^ -P ) 



(27r)"/' 

Note that we can write the brackets as (cf. eq. (20)) 



(48) 



CO 



2^0 (P - /i^ +iO)^ 



(49) 



For the spatial Fourier transform (48), we use Bochner theorem (cf. 
eq.l4)) 

^ f /„ 9\-il , s / si ©(— /' „ , "-1 sin \ t \ fl ^ , , , 



Formula 6.737-5, p.761 of the table of ref.[ll] gives (6 = i// + 0) 



And of course: 

T [{P ~ (x) ^ for xo>0 

The retarded Fourier transform reproduces (48), with a change of sign, 
for xo > and is zero for xo < 0. We then get, for the Wheeler Green 
function (46): 

22 



^{{P-f^V} = {p^Q'-) (51) 

Again, the Wheeler propagator has support inside the hght-cone. But, 
instead of a Bessel function of the first kind, we have now a Bessel function 
of the second kind. 

It is clear that by taking into account eq.(28), we can evaluate convolu- 
tions of different Green functions. As it was done in section 3 for bradyons 
fields. 

5 Fields with complex mass parameters 

The descomposition in Klein-Gordon factors of a higher order equation, of- 
ten leads to complex mass parameters. Equation (1) is an example. The 
constituent fields obey eq.(2). A simple higher order equation such as (3) 
presents the same behaviour. Of course for a real equation the masses come 
in complex conjugate pairs. We consider: 



(□ - M^) (f) = , M = m + ii^ (m > 0) (52) 

This type of equation has been analyzed elsewhere (ref.[5]). The Green 
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functions are inverses of P + —fP — pg, where 0=(p^ + M^)^/^. The two 
poles at po i move when varies from to oo, on a hne contained in a 
band of width ±i//, centered at the real axis. 

The retarded Green function is obtained with a po-integration that runs 
parallel to the real axis, with Impo > For the advanced solution, the 
integration runs below both poles {Impo < — |a*|): 



: / ji^-u ip.f ^^^Q( ^) ( int _ -int\ 





And according to ref.[ll] (6.735 - 5): 

For the retarded Green function we get the same answer, with 6(+i) 
replacing 9(— t). 

The Wheeler propagator is then: 
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j-|(p + m2) 'I (x) = |Mt-ig!^' ^Vi^ (mq!) (53) 

Now we have the general result: The Wheeler function propagates in- 
side the light-cone for any value of the mass, real (bradyons), imaginary 
(tachyons) or complex (M — ■m + i/i). 

In the case of complex masses, a natural definition for the Feynman prop- 
agator is obtained by a po-i^tegration along the real axis. It is not difficult 
to see that 

oo 


t sgnjj, — 1 J ^-3 irk) (54) 

The first term in the right hand side is the Wheeler function. The second 
term corresponds to a positive loop around the pole in the upper half-plane, 
and a negative loop around the pole in the lower half-plane. 

The conjugate Feynman function (not the complex conjugate one ), is 
obtained by changing the sign of the second term. I.e.: by changing the sign 
of both loops. 

The term in cosQ|t| can be read in ref.[ll] (6.735 - 6). With these defini- 
tions the Wheeler function is half Feynman and half its conjugate. 
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sinfl I t 



Note that for a real mass, /j, corresponds to the "small" negative imaginary 
part added to the mass in Feynman's definition. Accordingly, for real mass 
we take sgn/i — —1. 

Equation (54) can be utilized to define: 

oo /"> I I 

J-{5(P + Mm(a;) = -^./fr^ f dk k'-^'-^^J^irk) (55) 



(P + M^) ~^ ={p + M^) + i7r5 (P + M^) (56) 
The last formula is valid for any mass, real or complex. 

6 Associated vacuum 

It is well known that the perturbative solution to the quantum equation of 
motion leads to a Green function which is the vacuum expectation value of 
the chronological product of field operators (VEV). Furthermore, when the 
quanta are not allowed to have negative energies, the VEV turns out to be 
Feynman's propagator. 

However, when the energy-momentum vector is space-like the sign of its 
energy component is not Lorentz invariant. It is then natural to have symme- 
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try between positive and negative energies. It has been shown in references 
[7] and [8] that under this premise, the VEV is a Wheeler propagator. 

To see clearly the origin of the difference between both types of prop- 
agators, we are going to compare the usual procedure with the symmetric 
one. 

A quantized Klein-Gordon field can be written as: 



For simplicity, we are going to consider a single (discretized) degree of 
freedom. 

The raising and lowering operators obey: 




(57) 



where 




a, a 



.+ _ 



= 1 



(58) 



The energy operator is: 
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Usually, the energy is redefined to be ho. The vacuum then obeys: 



/io I >= (59) 
It is a consequence of (58) and (59) that: 

<0|aa+|0>=l , <0|a+a|0>=0 (60) 

On the other hand, the symmetric vacuum is defined to be the state that 
has zero "true energy" : 



/i|0 >= I (aa+ + a+a) |0 >= (61) 
Equations (58) and (61) imply: 



< 0|aa+|0 >= ^ , < 0|a+a|0 >= (62) 

Let as assume, for the sake of the argument, that we define a "ceiling" 
state (as opossed to a ground state): 

a+|0>=0 (63) 

Equations (58) and (63) give: 
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< 0|aa+|0 >= , < 0|a+a|0 >= -1 (64) 

The usual normal case, eq.(60) leads to the Feynman propagator. The 
"inverted" case, eq.(64), leads to its complex conjugate. Then eq.(62), which 
is one half of (60) and one half of (64), leads to one half of the Feynman 
function and one half of its conjugate. This is the Wheeler propagator defined 
in section 1. 

The space of states generated by succesive aphcations of a and a+ on the 
symmetric vacuum has an indefinite metric. 

The scalar product can be defined by means of the holomorphic represen- 
tation [14]. The functional space is formed by analytic functions f{z), with 
the scalar product: 

j dzdze-^^fiz)^ (65) 

Or, in polar coordinates: 

oo 2tt 

if, g)^ Jdp pe-"' j d<t>f{z)W) (66) 



The raising and lowering operators are represented by: 
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a'^ — z , a = 4- (67) 

dz 



The symmetric vacuum obeys: 

\dz dz J \ dz J 

whose normahzed solution is: 

The energy eigenfunctions are: 



fn 



-1/2 



z-^/^z" (68) 



7 Unitarity 

In QFT, the equations of motions for the states of a system of interacting 
fields are formally solved by means of the evolution operator. 

U{t,to) \to >= \t> 

The interactions between the quanta of the fields is supossed to take place 
in a limited region of space-time. The initial and final times can be taken to 
be to — )■ — oo and t — )■ +00. Thus defining the S-operator: 
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S = U{+oo, -oo) 

We do not intend to discuss the possible problems of such a definition. 
Here we are only interested in its relation to Wheeler propagators. 

Usually, the initial and final states are represented by free particles. How- 
ever, when Wheeler fields are present, their quanta either mediate interac- 
tions between other particles, or they end up at an absorber. This circun- 
stance had been pointed out by J.A.Wheeler and R.P.Feynman in references 
[1] and [3]. In consequence, the S-matrix not only contains the initial and 
final free particles, but it also contains the states of the absorbers. Through 
the latter we can determine the physical quantum numbers of the Wheeler 
virtual "asymptotic particles" . For these reasons, even if the initial and fi- 
nal states do not contain any Wheeler free particle, for the verification of 
perturbative unitarity it is necessary to take them into account. 

We shall ilustrate this point with some examples. Let us consider four 
scalar fields 0^ (s=l,...,4) obeying Klein-Gordon equations with mass pa- 
rameters ml and the interaction A = They can be written as in 
eq.(57). 
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Unitarity implies: 

SS+ = 1 

or, with S 

We introduce the initial and final states and also a complete descompo- 
sition of the unit operator: 

< a\T + r+l/3 >= j da^ < a\T\-f >< -f\T+\/3 > 
For the perturbative development: 

n 

n— 1 ,. 

< a\Tn + r+l/3 >= J2 da^< ^l^n-.lT >< > (69) 

s=l •' 

In particular, To = and Ti=pure imaginary. 
For n=2 

< a\T2 + T+\P >= J da^ < a\Ti\-f >< 7|T+|/3 > (70) 

where we will take Ti = i0i020304- 
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01 and 02 are supposed to be normal fields whose states can be obtained 
from the usual vacuum. 



\a 



>=aja{\0> , 1/5 >= aja+,|0 > 



On the other hand, for 03 and 04 we leave open the posibility of a choice 
between the usual vacuum or the symmetric one. 

The left hand side of (70) comes from the second order loop formed with 
the convolution of a propagator for 03 and another for 04. When both fields 
are normal, we have the convolution of two Feynman propagators, where the 
real part is: 



Re 



(P + - iO) % (P + - iO) ^ = (p + m^) % (P + mfj - 
n^S (P + m^) * 5 (P + mfj 



And according to (42), we have in the physical region (P negative): 



Re 



(p + ml- io) ^ * (P + ml-io) ^ = 2 (p + m^) ^ * (p + m^) ^ 



(P<0) 



(71) 



Equation (71) implies that the left hand side of (70) for Feynman particles 
is twice the value corresponding to Wheeler particles. 
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The relation (70) is known to be valid for normal fields. So, there is no 
point in proving each. We are going to show where the relative factor 2 comes 
from. 

The decomposition of unity for normal fields is: 

1 = j da^ |7 >< 7| = |0 >< 0| + j d''-^q;a+ (q) |0 >< 0\a{q) + 

J d^'-^qi d''^'q2^a+ (gl) a+ {q^) |0 >< 0|a (gl) a (gl) + .... (72) 
Then, for the Ti matrix we have: 

< q;|Ti|7 >=< 0|ai (p) 0i(a;)|O >< 0|a2 ) 4>2{x)\Q > x 

< 0\Mx)at (ga) |0 >< 0\Mx)at (q^) |0 > (73) 
where an integration over x-space is understood. 

When the fields are expressed in terms of the operators a{q) and a'^{q), 
as in equation (59) we obtain: 



{2nr 6{p-qs- 



And 



{2nY S(^qs + q4-p') 

lr)\2{v-l) . I t t 
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Multiplying together (74) and (75) and adding all possible I7 >< 7I (all 
q2 and q^), we get: 



da^ < a\Ti\j >< j\T+\P > 
dq 



5(p- p') 



16 (27r)^^ ^Ju}iU}2U}'iU}'2 
S -UJs {(f) - oj.i{p- q)) 



(76) 



UJ3{q)uj4{p- q) 

This result coincides with (70) (l.h.s.) when the p°-convolution is carried 
out. 

Suposse now that one of the fields, says 04, has the Wheeler function as 
propagator. Instead of eq.(71) we have: 



Re 



[P + ml- iOj ^ * (^P + ml-io 

(^p + miy\ (^p + miy^ 



-1 



(77) 



Half the value of (71). 

To evaluate the matrix < Ti > for this case we note that the descompo- 
sition of unity for the states of ^4 (with an indefinite metric) is now: 

I = J day I7 >< 7I = |0 >< 0| + J d^'-^q V2a+ (q) |0 >< 0|a (g) V2- 



j d^-^q V2a{q) |0 >< 0|a+ {q)V2+ .... 



(78) 
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The normalization factors come from the VEV quoted in section 6, eq.(62). 
It is not necessary to evaluate again the matrix element (73). Its last vacuum 
expectation value has now a factor 1/2 from eq.(62), and a factor -\/2 form 
the normalization in (78). When the matrix for Ti and are multiplied 
together, we get an extra factor (-\/2/2)^= 1/2 As it should be for unitarity 
to hold. 

When both fields 03 and 04, are of the Wheeler type, we get for the 
convolution of the respective Wheeler propagators the same result (77). 

The matrix element of Ti gains now two factors \/2/2, i.e. a factor 1/2. 
When we multiply < Ti >< > we get a factor 1/2 • 1/2= 1/4. And we 
seem to be in trouble with unitarity. However, in this case a new matrix 
contributes to < Ti >. It is: 

< 0|ai (pi) (t)i{x)al {qi) [q'^ |0 >< 0|a2 (P2) 02(a;)a^ (52) ^2 (^2) |0 > 

< O|03(a;)a+ (gs) |0 >< O|04(a;)a+ (^4) |0 > (79) 

(79) is only possible when both, ^3 and 04 are associated with symmetric 
vacua. 

For the first matrix factor wc have: 

< 0|ai (pi) (pi{x)a]; {qi) a]; (q^ |0 >= 
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A similar matrix factor from < > gives: 



< 0|ai (qi) ai (g^) 0i(y)ai (p^) |0 >= 



S (Pi - ?i) , — + S(p,- gi) , ^ ^ (81) 



When we multiply together (80) and (81), the crossed terms do not con- 
tribute {S{pi — p'l) — 0). The other two terms give equal contributions. A 
similar evaluation can be done for the second factor of (79) and the corre- 
sponding factor of < T+ >. For this reason we are going to keep only the 
first terms from (80) and (81) (multiplied with the appropriate constants): 



2 Q-ili^ 

<a\Ti\^>= 5{p,-q,) 5{p2-q2) x 

(27r) ^2uji (gi) 



^2a;2 (gs) V2^3 (gs) 2^2a;4 (^4) 
And after performing the x-integration, 

,^ , (27rY S (-Qi - ?2 + 93 + qi) S (pi - qi) 6 (pa - ^2) 

< Q;|ri|7 >= 2(l/-i) , / ^ 

2(27r) A^Juj[uj'2UJzuJa 

Analogously: 
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{2nY HQi + Q2-q3-q4)S(p[-q[)s(p'2-q'2) 
< l\Ti 1/3 >= 2(.-i) I , , 

The sum / da-y < a\Ti\j >< j\Ti~\(3 > corresponds to an integration on 
(fi,(f{,(f2,92- is see that after this operations we get one fourth of 

(76). Thus completing the proof of unitarity, for the proposed example. 

Let us now consider the case in which 03 and 04 are fields obeying Klein- 
Gordon equations with complex mass parameters (See section 5). 

The solution of eq.(52): 



(□ - M^) = 0, M ^m + in (m > 0) 



can be written as: 



(27r)~ J 



p 



a {p) e-'^* + b{p)e 



iQt 



where 



The quantization of leads to the rules (ref.[5]): 



(82) 



(83) 



{P),b{p')] =S (p-p) 
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and to the adoption of the symmetric vacuum. From which we get: 



< 0\a ip) b {p') |0 >= - < 0|6 (p) a {p') |0 >= ^5 (p- p) (84) 

Accordingly, the descomposition of unity for complex mass fields is (com- 
pare with (78)): 

I = |0 >< 0| + j (T'^q V2h (q) |0 >< 0|a (q) |0 > \/2- 

j (f-^q \/2a (q) |0 >< 0|6 {q)\0>V2 + .... 

The matrix < q;|Ti|7 > have the form given in eq.(73), but now, instead 
of 03" and 04 , we have to write the four possible operators: 03, 04; 03, 64; 63, 04 
and 63, 64. 

When multiphed with < ^\T^\/3 > as in (76) they give similar contribu- 
tions except for the signs of 003 and 004 in the arguments of the 5-functions. 
Then, the 5-function for each of the terms in the integral of (76) (r.h.s), should 
he:S{po — 0,3 — 0,4), S{po — ^3 + 0,4), ^{Po + ^^3 - ^4) and S{po + ^3 + ^4): 
respectively. 

This is exactly one half of the convolution product of the Wheeler prop- 
agators for 03 and 04 . The other half comes from the contribution of the 
matrix elements of the form (79). 
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Of course, for every field corresponding to a complex mass M, there is 
another field corresponding to the complex conjugate mass M*. 
We have: 



< 0|6+ (p) a+ (p) |0 >= - < 0|a+ (p) 6+ (p') |0 >= ^5 (p -p' 
I = |0 >< 0| + y" d'^-^g y2a+ {q) |0 >< 0|6+ (g) |0 > ^/2- 
j d^-^q \f2h-^ {q) |0 >< 0|a+ {q) |0 > .... 

The T-matrix can be constructed with the four possible contributions: 
03,04; 03, 01; 03 , 04 and 0^, (pt- 

Correspondingly, for each Wheeler propagator with mass M, there is an- 
other one with mass M*. And of course, four possible convolution products. 
It is easy to see that the total convolution is real, as well as the total T- 
matrix. 

The proof of unitarity for the chosen example, where the masses are com- 
plex numbers, contains as particular cases all fields with symmetric vacua. 
For bradyons M=m. For tachyons, the hmit M — > i/i, must be taken. 
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Similarly in other cases a proof of unitarity for Feynman propagators, 
based on the decomposition (72), can be converted into a proof of unitarity 
for Wheeler propagators, by using the corresponding decomposition (78). 

8 Discussion 

We have shown that the Wheeler propagator has several interesting proper- 
ties. In the first place we have the fact that it implies only virtual propa- 
gation. The on-shell 5-f unction, solution of the free equation, is absent. No 
quantum of the field can be found in a free state. The function is allways 
zero for space-hke distances. The field propagation takes place inside the 
light-cone. This is true for bradyons, but is also true for fields that obey 
Klein-Gordon equations with the wrong sign of the mass term and even for 
complex mass fields. The convolution of two Wheeler functions gives an- 
other Wheeler function. In p-space, this convolution coincides for P < 0, 
with the convolution of the two on-shell (5-functions. In spite of the fact that 
each of the latter only contains free propagation, while each of the former 
only contains virtual propagation. The usual vacuum state is annihilated 
by the descending operator a, and gives rise to the Feynman propagator. 
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The Wheeler Green function is associated to the symmetric vacuum. This 
vacuum is not annihilated by a, but rather by the "true energy" operator, 
a symmetric combination of annihilation and creation operators. The space 
of states generated by a and a"*" has an indefinite metric. There are known 
methods to deal with this kind of space. In particular we can define and 
handle all scalar products by means of the "holomorphic representation" [14] . 
Due to the absence of asymptotic free waves, no Wheeler particle will appear 
in external legs of the Feynman diagrams. Only the propagator will appear 
explicitly, associated to internal lines. So the theoretical tools to deal with 
matrix elements in spaces with indefinite metric, will not, in actual facts be 
necessary for the evaluation of cross-sections. However, the descomposition 
of unity for spaces with indefinite metric, is needed for the proof of another 
important point. The inclusion of Wheeler fields and the correspondings 
Wheeler propagators do not produce any violation of unitarity, when only 
normal particles are found in external legs of Feynman diagrams. 

To complete the theoretical framework for a rigorous mathematical anal- 
ysis, it is perhaps convenient to notice that the propagators we have defined, 
are continuous linear functionals on the space of the entire analytic functions 
rapidly decreasing on the real axis. They are known in general as "Tempered 
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Ultradistributions" [15, 16, 17, 18]. The Fourier transformed space contains 
the usual distributions and also admits exponentially increasing functions 
(distributions of the exponential type) (see also ref.[19]). 

We must also answer the important question. What are the possible uses 
of the Wheeler propagators?. 

In the first place we would like to stress the fact that the quanta of 
Wheeler fields can not appear as free particles. They can only be detected 
as virtual mediators of interactions. It is in the light of this observation that 
we must look for probable apphcations. 

We will first take the case of a tachyon field. It is known that unitar- 
ity can not be achieved, provided we accept the implicit premise that only 
Feynman propagators are to be used, with the consequent presence of free 
tachyons. This work can also be considered to be a proof of the incompatibil- 
ity of unitarity and Feynman propagator for tachyons. To this observation 
we add the fact that if the propagator is a Wheeler function, a tachyon 
can not propagate freely. Consequently we are led to the acceptance of the 
complete spectrum: < fj,"^ and > /x^. With the remark that the real 
exponentials one gets for p ^ < //^ are not eigenfunctions of the Hamiltonian 
(ref.[8, 20]). Furthermore, this procedure fits naturally into the treatment 
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for complex mass fields of section 5. To this case it could be related the 
Higgs problems. The scalar field of the standard model behaves as a tachyon 
field for low amplitudes. The fact that the Higgs has not yet been observed 
suggest the posibiUty that the corresponding propagator might be a Wheeler 
function [21]. It is easy to see that this assumption does not spoil any of the 
experimental confirmations of the model (On the contrary, it adds the non 
observation of the free Higgs) . 

Another possible application appear in higher order equations. Those 
equations appear for example in some supersymmetric models for higher di- 
mensional spaces [6]. They can be decompossed into Klein-Gordon factors 
with general mass parmeters. The corresponding fields have Wheeler func- 
tions as propagators. It is interesting that there are models of higher order 
equations, coupled to electromagnetism, which can be shown to be unitary, 
no matter how hig the order is [22]. 

The acceptance of thachyons as Wheeler particles, might be of interest for 
the bosonic string theory. With the symmetric vacuum we can show that the 
Virasoro algebra turns out to be free of anomalies in spaces of arbitrary num- 
ber of dimensions [23]. The excitations of the string are Wheeler functions 
in this case. 
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